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Abstract 

The search of the correct equation of motion for a classical charged particle under the action 
of its electromagnetic (EM) self-field, the so-called radiation-reaction equation of motion, remains 
elusive to date. In this paper we intend to point out why this is so. The discussion is based on 
the direct construction of the EM self-potentials produced by a charged spherical particle under 
the action of an external EM force. In particular we intend to analyze basic features of the LAD 
(Lorentz-Abraham-Dirac) and the LL (Landau-Lifschitz) equations. Both are shown to lead to 
incorrect or incomplete results. 
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I. INTRODUCTION 



An ubiquitous phenomenon which characterizes the dynamics of classical charged par- 
ticles is the occurrence of self-forces, in particular the electromagnetic (EM) one which is 
produced by the EM fields generated by the same particles. It is well-known that such a 
force acts on a charged particle when it is subject also to the action of an arbitrary external 
force (see Lorentz, 1892 see also for example Landau and Lifschitz, 1957 [2I). This 
phenomenon is usually called as radiation reaction (RR) (Pauli [4]) or radiation damping 
(see (3]), although a proper distinction between the two terms can actually be made [fj. In 
the sequel we denote as RR problem the treatment of the classical mechanics of a charged 
particle in the presence of its EM self-field, and RR equation the corresponding equation of 
motion for such a particle. The RR problem, apart its almost endless physical applications, 
is of fundamental importance from the theoretical viewpoint for the formulation of consistent 
relativistic theories, in particular kinetic theory, plasma dynamics and astrophysics. Despite 
efforts more than a century-long, the problem of its " exact" theoretical description remains 
still elusive. This refers, in particular, to the construction of a relativistic equation of motion 
which results both non-perturbative, in the sense that it does not rely on a perturbative ex- 
pansion for the electromagnetic field generated by the charged particle and non- asymptotic, 
i.e., it does not depend on any infinitesimal parameter (such, for example, as the radius 
of the particle, if it is identified with a spherically-symmetric charge distribution of finite 
radius, as in the Lorentz approach |l|). The purpose of this paper is to point out that such a 
problem, since the historical work of Lorentz which first investigated it [1], still remains es- 
sentially unsolved to date. To prove the above statement, in this paper we intend to analyze 
in detail the derivation and basic features of the relativistic equations of motion available 
so-far, with particular reference to the so-called LAD and to the LL equations, respectively 
named after Lorentz, Abraham and Dirac and Landau and Lifschitz |2J. A side motivation 
of this work is also provided by the recent claim (Rohrlich, 2001 see also Spohn, 2000 
[3]) that the LL equation should be regarded as the exact (i.e., both non-perturbative and 
non-asymptotic) RR equation. 
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II. A RE-FORMULATION OF THE TRADITIONAL APPROACH TO THE RR 
PROBLEM 



In the weakly-relativistic treatment of classical mechanics, the traditional form of the 
or a charged particle subject to its own EM self-field can be obtained 



equations of motion 

(see for example Ref. 2|) by a suitable asymptotic expansion of the EM self-potentials which 
generate such a force. It is instructive to discuss here the basic steps of the derivation. 
For definiteness, we shall consider here the case of a single point-particle with distributed 
charge density (here denoted as finite-size charge) carrying constant mass and charge, m 
and q (this choice is actually analogous to that adopted in the original Lorentz approach 
More precisely, it is assumed that the particle mass is concentrated in the center of the 
sphere (so that the particle degree of freedom is the same as that of a point particle), while 
its charge is assumed to be uniformly distributed on a spherical shell of finite radius o > 0, 
which carries the homogeneous surface charge density p = q/ATra 2 (finite-size spherical- 
shell charge). Thus, letting R = |r — r'| and denoting respectively [p(v', t), J(r', t)] and 
\p(r',t — J(r', t — ^)] the charge and current densities of the particle at time t and at 
the retarded time if = t — — , we shall impose that 

P(r,t) = -^5(\r-r(t)\-a), (1) 
J(r,f) = ^- 2 r(t)5(\r-r(t)\-a). (2) 

Then, if r(t) is at time t the position of the particle center of mass, the EM self-force acting 
on the same particle is represented by the Lorentz force acting on the particle itself, i.e., 



F {self \r(t) } r(t) } t)= J d 3 r' 



p(r',t)E( seZ/ V,t) + ^J(r',t) x B< aeV V,t) 



(3) 



where F( se '^ is usually denoted as RR reaction force acting on the particle. Here 
{E* 56 ^, B( se/ ^} denotes the EM self-field generated by the particle itself, related to 
corresponding EM potentials |0( se/ /) ; A^ se/ ^^} , here denoted as EM self-potentials. By 
definition |0( se/ /) ; A^ se/ ^} - to be evaluated at time t and at the position r [to be 
later identified with the particle center of mass r(t)] - can be defined as 0^ se/ ^^(r, t) = 
Jd 3 r'^ [p(r',t- f) -p(r',t)] and A^{r,t) = \ J d 3 r'^ [J(r', t - f ) - J(r', t)] . We 
stress that the second terms on the r.h.s. of these equations take into account the charge 
and current densities evaluated at the same time t (current time), which manifestly do not 



contribute to the self potentials when they are evaluated at the particle position. In order 
to evaluate explicitly the EM self-potentials in this case [Eqs.flT]), (j2J)], let us now introduce 
for them an asymptotic expansion in power series of (3 = v/c, assuming that v/c <C 1. 
Retaining only terms up to third order (in 0) in (p^ sel ^(r,t) and respectively first order (in 
(5) in A( sel f'(r,t), there results 



cdt 

A^(r,t) = -~ I tfV'J(r',t). 



Ll(self) 



6c 3 dt 3 



(5) 



self) + -§-J = 0, A'( se// )(r,t) = A( se// )(r,t)-V/, 



Thus, introducing the EM gauge r — Y 
where / is the EM gauge function (to be denoted for later reference as RR gauge) 



(6) 



the transformed vector self-potential A'( sel f\r,t) becomes 

A^(r,t)^-~ j dVj(r',t)+V J d 3 r'p( r >,t)-~V J d 3 r'Rp(r',t)+-^^V j d 3 r'R 2 p(r',t). 

(7) 

After straightforward calculations one obtains 

f r(t) 



F^\r(t) } r(t),t) = -^r(t) + aq 2 



+ 



era 



2a 3 



r{t)/3 2 



(8) 



where f3 2 = r 2 (t)/c 2 and in the case of the finite-size spherical shell here considered there 



results a = 1. Thus, one recovers in this way t 



re weakly-relativistic RR equation or the 



weakly-relativistic LAD equation (after Lorentz {2], Abraham 8] and Dirac [9|) of the form 



IJIrY = F + g. 



(9) 



Here the notation is standard. Thus, 



2q 2 - aq 2 2 
r + -^—^r{t)p 



3c 3 2a 3 

is the weakly-relativistic EM self-force. Moreover 



q 



g(art) _|_ l r x -Q(ext) 



(10) 



m R = m- m EM , 



m EM = az-^, 



(11) 



are respectively the Lorentz force produced by the external EM field, the renormalized mass 
and the so-called EM mass. Eq.^j manifestly indicates that the limit a —>■ (point-particle) 
does not exist. Hence the LAD equation is only defined for a finite-size spherical- shell charge, 
i.e., for which it results o > and possibly also > (the vanishing of involves in fact 
the manifest violation of Newton's second law) . The derivation of the analogous relativistic 
RR equation is usually achieved (see for example the treatment given in Ref. [21]) under the 
implicit assumptions that: 1) higher order corrections in (3 are negligible; 2) all previous 
asymptotic expansions in power series of (3 remain valid for the whole range of values of 
particle velocities, i.e., even if (3 ~ o(l), namely arbitrary close to the speed of light. In 
practice, the construction of the equation is achieved formally by writing Eq. in covariant 
form and by suitably replacing the definitions of the renormalized mass m# and of the self- 
force g by an appropriate Lorentz scalar m R and a 4- vector g^. Both are to be defined 
in such a way that in the weakly-relativistic limit they recover the correct values set by 
the previous equations. The relativistic RR equation of motion thus obtained, nowadays 
popularly known as the relativistic LAD equation, 

m oR c— = q -F^u v + g» (12) 
as* c 

was first presented by Dirac in his famous paper on relativistic radiation reaction in classical 
electrodynamics jjj. Here, g^ is the EM self- force which reads 



Furthermore the relativistic renormalized mass has the form 

m R = m + m EM + Am EM , (14) 

with m , rriEM an d Attiem respectively , the rest mass, the EM mass and a suitable relativis- 
tic correction. Finally, as usual r M , and F^ u denote the 4-position vector, the 4-velocity 
vector and the Faraday tensor associated to the external EM field A^ ext '^. Nevertheless, 
the LAD equation is valid only in the case of the Minkowsky metric. Its gene ralization for 



arbitrary curved space-time was later carried out by DeWitt and Brehme 10] 
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III. DIFFICULTIES WITH RR EQUATIONS 

Since Lorentz famous paper [l| many textbooks and research articles have appeared on 



the subject of RR. Among them are 



HQ, 



14 
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161 ]. where one can find the discussion 



of the related problems: mass renormalization, non-uniqueness, runaway solutions and pre- 
acceleration. In the literature the usual derivations of relativistic form of the EM self-force 
are made under the implicit assumption that all the expansions in powers used near the 
particle trajectory are valid for the whole range of values of particle velocity, in particular, 
arbitrary close to that of the light. But it is easy to see that this is not true in general. 

It is often said that the LAD equation, in its weakly- relativistic form given by Eq.([5]), is 
unsatisfactory because it requires the specification of the initial acceleration r(t D ), besides 
the initial state |r(t D ), r(t Q )| . As a consequence it violates the Newton's principle of deter- 
minacy (NPD), one of the building blocks of classical mechanics. Moreover, another critical 
issue arises due to the appearance of so-called runaway solutions. These are solutions 
which by definitions blow up in time for t — > +oo. One can show that they are present for 
arbitrary external forces F. In particular, in the case F = the general solution of Eq.([n]) 
takes the form r(t)=r(t ) exp | ft~ to ) j , where r > is the constant parameter r = J^s- This 
equation is clearly unsatisfactory. In fact, it violates also another fundamental principle of 
classical mechanics, the Galilei law of inertia, according to which an isolated particle must 
have a constant velocity in any inertial Galilean frame. In an attempt to circumvent this 
difficulty it has been suggested to replace Eq.Q by a suitably-modified integral (or differ- 
ential) equation. An example is provided by an integro-differential equation which in its 
weakly-relativistic form reads 12 , 



m R r = J F(t + rs)e~ s ds = -J F(x)e~^dx (15) 
o t 

(Eaag equation). This equation is manifestly consistent with PND and the law of inertia 
(since runaway solutions are excluded). However, it exhibits another serious drawback: it 
violates the principle of causality. Indeed, contrary to the requirement according to which 
the effect cannot precede the cause, the paradox of pre- acceleration is implied by Eq.f U5l) . 
Accordingly, a particle obeying Eq. ffT5l) would experience a force before actually turning it 
on! Eq. (fT5|) clearly violates this principle since if a constant force F Q is suddenly turned on 
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at t D letting F(x) = for x < t Q and F(x) = F G for x > t Q (sudden force) it follows 

oo 

'dx = F G - / e~~dx = F a e~ r t < t a , 

(16) 

dx = F a t > t a . 

Hence, the RR equation before and after t Q will be respectively provided by 

, F e~^ t < t , 
m R r = \ (17) 

F t > t a . 

Hence, the particle should actually experience an acceleration for all t < t Q , a phenomenon 
usually denoted as pre-acceleration. For this reason neither Eq.Q nor (fT5]) can be considered 
satisfactory RR equations. 

Another attempt is to adopt an "iterative approach" whereby the time derivative of the 
acceleration appearing in the self-force (flOj) is approximated in terms of the time derivative 
of the external force (this approach is manifestly admissible only if the self-force can be 
treated as suitably small). In the weakly-relativistic LAD, letting first ^§-r(t)r = t^F, 
one obtains the reduced equation m^r = F+r^F. More generally, by carrying out the full 

oo 

iteration in Eq.([9]), one obtains the equation m R r = F+ r n ^-F. This equation, and the 

n=l 

analogous relativistic one which follows in a similar way from Eq. (fl2|) . are however useless 
since they are infinite-order ode's (and hence have no solutions!). 

Nevertheless, this difficult can be circumvented by introducing a "reduction method" 
whereby all the derivatives higher than r appearing in each term 4^F are expressed via 
the same iterative approach which permits to cast -CrF. for each n > 1 in the form 



j|_F = F n (r, r, t). This idea has been exploited by Cook [18( who obtained the so-called 



Cook RR equation 

oo 

m H r = F+^r"F n (r,r,t). (18) 
This procedure was actually first adopted by Landau and Lifschitz |2J based on a one-step 
iteration of the form 

m R r ^ F+rFi(r,r,t) (19) 

(which could be denoted as the weakly-relativistic LL equation |2j). In the corresponding rel- 
ativistic formulation (see Ref. [2j), applying the one-step reduction process to the relativistic 
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LAD equation (|T2|) delivers for [as given by Eq. (fl3l) ] the approximation 



2m R & 



Q in- 



The relativistic LL equation obtained in such a way by Landau and Lifschitz 
volves, however, also the replacement of the renormalized mass m R with the inertial mass 
only m . While in the framework on classical electrodynamics the latter position remains 
totally unjustified, the resulting equation has been claimed by Rohrlich [6j) to be the ex- 
act RR equation. Also in view of the previous discussion, the validity of this statement 
seems unlikely. In fact the LL equation (both in the weakly-relativistic and fully relativistic 
versions) has remaining serious problems: 

• One reason is that as recalled in Sec. 2 the theory is asymptotic in (5 and does not take 
into account in a consistent way higher-order finite-/? effects; 

• A second motivation is that it only applies provided the external force is suitably 
smooth, i.e., of class C 2 . In fact, if an external force is turned on suddenly (sudden 
force), as in the example provided above for the weakly-relativistic case, the LL equa- 
tion becomes manifestly invalid. Since sudden forces cannot be ruled out, purely on 
first principles, this is actually a major conceptual difficulty (and potential inconsis- 
tency) of the part of the LL equation. 

• An additional difficulty of the LL equation is that it is non-variational, i.e., the Hamil- 
ton variational principle does not apply for such an equation, even if the external force 
is identified with the Lorentz force [see Eq. lTTTl) ]. This feature is actually common to 
all current RR equations. This is illustrated, for example, by the weakly-relativistic 
LAD equation Eq.Q which appears manifestly non-variational. In fact in this case the 
Lagrangian function L should actually depend not only on the Lagrangian state (r,r) 

but also on the acceleration r. Thus, for example, letting introducing the Lagrangian 

2 r i 
L(r,r, r, t) = \mRY — q <p^ ext \r,t) — -r • A^ ez ^(r,t) + /(r,r, r, t), for the validity of 



S 



the Hamilton variational principle it should be 

*2 *2 



5 J dtL(r(t),r(t),'r(t),t) = J dtSr(t) ■ m R r - f] + (21) 

ti 

+ y dt6f(r(t)j(t),i(t),t) = 0, 



ti 



where 8f(r(t),r(t),r(t),t) is an exact differential form such that 

ti ti 2 

/ dtSf(r(t),r(t),r(t),t) = J dtSr(t) • f^r(t). This implies, however, 

h ti 



.2 h 



j dt8f(r(t),r(t), r(t),t) = J dt [8r(t) ■ r(t) - 5r(t) ■ r(t)] = 0, (22) 

which means that a real function f(r(t),r(t),r(t),t) fulfilling the previous constraint 
cannot exist. 

• Finally, another serious difficulty is related also to the conditions validity of the reduc- 
tion process indicated above. In fact, it is obvious that the one-step reduction adopted 
for the derivation of the LL equation provides, at most, an asymptotic approximation 
for the (still unknown) exact RR equation. 



IV. CONCLUSIONS 



In this paper we have reviewed aspects of the RR problem. We have shown that: 

• the difficulties with the LAD equation are intrinsic, i.e., arise due to the fact that the 
equation is a third order ode. As a consequence, the Newton's principle of determinacy 
and the Galilei law of inertia are potentially violated; 

• attempts to circumvent these difficulties, based on various approaches (Haag, LL and 
Cook equations) fail for different reasons; 

• limitations of the current RR equations have been pointed out, with particular refer- 
ence to the conditions of validity of the LL equation; 
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all current RR equations are non-variational, implying - contrary to common knowl- 
edge in classical mechanics - that the dynamics of a charged particle described by 
these model equations does not define an Hamiltonian system. This is the reason 
why relativistic systems of charged particle usually are not (or cannot be) described 
as Hamiltonian systems. Nevertheless, it not clear whether this feature is only an 
accident, i.e., is only due to the approximations introduced in the RR equation or is 
actually intrinsic to the nature of the RR problem. Unfortunately, the precise solution 
of this dilemma is not yet known, although the prevailing opinion seems directed to 



the first possibility (see for example Ref. [19(]). For a proper discussion of the issue we 
refer to the accompanying paper j^ . 

This suggests, in our view, that in many respects the RR problem is still open. Its so- 
lution should be based on the search of a relativistic, non-perturbative equation of motion 
for a particle in the presence of its EM self-field. This problem poses not just an intellec- 
tual challenge but a fundamental physical requirement in all applications which involve the 
description of relativistic dynamics of classical charged particles, such as relativistic kinetic 
theory, plasma physics and astrophysics. In an accompanying paper an attempt at a possible 



solution of the RR problem will be discussed 



20|. 
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